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ABSTRACT 

We present a self-consistent, absolute isochronal age scale for young (< 200Myr), 
nearby (< 100 pc) moving groups in the solar neighbourhood based on homogeneous 
fitting of semi-empirical pre-main-sequence model isochrones using the r 2 maximum- 
likelihood fitting statistic of |Naylor fe Jeffries] in the My,V — J colour-magnitude 
diagram. The final adopted ages for the groups are: 149ljg Myr for the AB Dor moving 
group, 24 ± 3 Myr for the Pic moving group (BPMG), 45^7 Myr for the Carina 
association, 421® Myr for the Columba association, 11 ± 3 Myr for the r/ Cha cluster, 
45±4Myr for the Tucana-Horologium moving group (Tuc-Hor), 10±3Myr for the TW 
Hya association, and 22lg Myr for the 32 Ori group. At this stage we are uncomfortable 
assigning a final, unambiguous age to the Argus association as our membership list 
for the association appears to suffer from a high level of contamination, and therefore 
it remains unclear whether these stars represent a single population of coeval stars. 

Our isochronal ages for both the BPMG and Tuc-Hor are consistent with recent 
lithium depletion boundary (LDB) ages, which unlike isochronal ages, are relatively 
insensitive to the choice of low-mass evolutionary models. This consistency between 
the isochronal and LDB ages instills confidence that our self-consistent, absolute age 
scale for young, nearby moving groups is robust, and hence we suggest that these ages 
be adopted for future studies of these groups. 

Software implementing the methods described in this study is available from http: 
//www.astro.ex.ac. uk/ people/timn/tau-squared/. 

Key words: stars: evolution - stars: formation - stars: pre-main-sequence - stars: 
fundamental parameters - techniques: photometric - solar neighbourhood - open clus¬ 
ters and associations: general Hertzsprung-Russell and colour-magnitude diagrams 


1 INTRODUCTION 


Over the past couple of decades, hundreds of young low- and 
intermediate-mass stars have been discovered in close prox¬ 
imity to the Sun. These stars are not uniformly dispersed 
across the sky, but instead comprise sparse, (mostly) gravita¬ 
tionally unbound stellar associations within which the mem¬ 
bers share a common space motion. Approximately 10 such 
moving groups, with ages of between ~ 10 and 200 Myr, have 
been identified within a distance of 100 pc from the Sun (see 
e.g. |Zuckerman &; Song||2004| |Torres et al.||2008] |Mamaje~k| 


2015). Given their proximity to Earth, the members of these 


groups therefore play a crucial role in our understanding of 
the early evolution of low- and intermediate-mass stars. Fur¬ 
thermore, such stars provide ideal targets for direct imag- 


* E-mail: cbell@pas.rochester.edu (CPMB) 


ing and other measurements of dusty debris discs, substellar 


objects and, of course, extrasolar planets (e.g. Janson et al. 


2013 Chauvin et al.|2015 MacGregor et al.|2015l. 


Whilst the relative ages of these young associations are 
well known e.g. the TW Hya association (TWA) is younger 
than the /3 Pic moving group (BPMG), which in turn is 
younger than the AB Dor moving group, the absolute ages 
of these groups are still under-constrained. There are sev¬ 
eral methods for deriving the absolute age of a given young 
group of stars (see the review of Soderbl om et alT| |2014| ) , 
however there are known problems with almost every one. 
Arguably the most common age-dating technique, using the¬ 
oretical model isochrones, still suffers from a high level of 
model dependency which arises from the differences in the 
treatment of various physical aspects, as well as the val¬ 
ues of adopted parameters; most importantly the treatment 
of convection, the sources of opacity, the handling of the 
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stellar interior/atmospheric boundary conditions and the 


initial chemical composition (see e.g. Dahm 2005 Hillen- 


brand, Bauermeister, & White|2008l. Even model indepen¬ 


dent methods, such as using kinematic information to in¬ 
fer the expansion rate of a given group or estimating the 
time at which said group occupied the smallest volume in 
space, suffer from issues of subjectivity with regard to the 


stars which are included/excluded (e.g. Song, Zuckerman, & 


Besseh|2003j |de la Reza, Jilinski, fe Ortega|2 006| lSoderblom 

et al.|2014| ). Furthermore, these kinematic methods have also 
been shown to be unreproducible when similar analyses have 


been performed using improved astrometric data (e.g. Mur- 


phy, Lawson, fe Bessell||2013 Mamajek & Bell 2014). 


Despite the apparent problems with deriving absolute 
ages for young associations, recent age estimates based 
on the lithium depletion boundary (LDB) technique have 
started to instill confidence that we have both precise and 
accurate ages (to within just a few Myr) for at least two of 
these groups - the BPMG and Tucana-Horologium moving 
group (Tuc-Hor). The LDB technique works by identifying 
the lowest luminosity within a given (presumed coeval) pop¬ 
ulation of stars at which the resonant lithium (Li) feature at 
6708 A shows that Li remains unburned. Although the LDB 
technique relies on the same theoretical models of stellar evo¬ 
lution from which (model-dependent) isochrones are created, 
the luminosity at which the transition from burned to un¬ 
burned Li occurs is remarkably insensitive to the inputs and 


assumptions noted above (see e.g. Burke, Pinsonneault, & 


Sills 2004 Tognelli, Prada Moroni, & D egrinnocenti|20T5 1. 
Given this high level of model insensitivity, [Soderblom| 


et al. (2014) propose that LDB ages therefore provide the 


best means of establishing a reliable and robust age scale for 
young (< 200 Myr) stellar populations. Unfortunately, there 
is a lower age limit to the applicability of the LDB tech¬ 
nique (~ 20 Myr) which arises due to a much higher level of 
model dependency between the various evolutionary mod¬ 
els in this age regime. Furthermore, for many of the young 
groups studied here, the censuses of low-mass members is far 
from complete and therefore calculating a robust LDB age 
is not possible. Therefore, if we are to establish an absolute 
age scale for young moving groups in the solar neighbour¬ 
hood, we must use different age-dating techniques, however 
we must also ensure that the resultant ages are consistent 
with those calculated using the LDB technique. 

In a series of papers ( |Bell et al.|2012||2013|[20l4) ) we dis¬ 
cussed some of the main issues with using pre-main-sequence 
(pre-MS) model isochrones to derive ages for young clusters 
in colour-magnitude diagrams (CMDs), and in particular, 
the inability of these models to reproduce the observed loci 
in CMD space for clusters with well-constrained ages and 
distances. In these papers we introduced a method of cre¬ 
ating semi-empirical pre-MS model isochrones using the ob¬ 
served colours of stars in the Pleiades to derive empirical 
corrections to the theoretical colour-effective temperature 
(T e tf) relations and bolometric corrections (BCs; hereafter 
referred to together as BC-T e ff relations) predicted by at¬ 
mospheric models. Of the clusters studied in the Bell ct al. 
studies, Jeffries et al. (2013) recently identified the LDB in 


NGC 1960 and derived an age of 22 ± 4 Myr (cf. with the 
isochronal age of 20 ± lMyr). The consistency between the 
isochronal and LDB ages lends confidence in the use of such 


semi-empirical model isochrones to establish an absolute age 
scale for other young clusters and associations. 

In this study we use four different sets of semi-empirical 
pre-MS model isochrones (which also include the effects of 
binary stars) in conjunction with a maximum-likelihood fit¬ 
ting statistic to derive absolute ages for eight young moving 
groups within 100 pc of the Sun. The ages of these groups 
are on the same scale as those derived for a set of young 


galactic clusters in Bell et al. (20131, and thus we are begin¬ 
ning to establish a robust set of ages for young benchmark 
stellar populations out to distances of ~ lkpc. In Section [2] 
we discuss our sample of young groups, and having compiled 
a list of members and high-probability candidate members 
for each, then collect broadband photometry and assign dis¬ 
tances. Section [3] introduces the theoretical stellar interior 
and atmospheric models, as well as providing a brief descrip¬ 
tion of our method for creating the semi-empirical model 
isochrones. In Section[4]we describe the maximum-likelihood 
fitting statistic and derive best-fit ages for each of the young 
moving groups in our sample. We then discuss our results in 
Section [5] and finally present our conclusions in Section [6] 


2 THE DATA 

We focus on the following young (< 200Myr), nearby (< 
100 pc) moving groups: the AB Dor moving group, Argus 
association, BPMG, Carina association, Columba associa¬ 
tion, r/ Cha cluster, Tuc-Hor, TWA, and 32 Ori group. As¬ 
signing members to young moving groups is typically based 
on a combination of kinematic diagnostics (e.g. proper mo¬ 
tions, radial velocities, etc.) in conjunction with youth indi¬ 
cators such as Ha and X-ray emission, and Li absorption. 
Having established a list of members and liigh-probability 
candidate members (hereafter simply referred to together as 
‘members’) for each group, we then compile broadband pho¬ 
tometric measurements and assign distances to these stars 
so that we can then place them in a CMD before fitting 
with model isochrones. Note that in our adopted member¬ 
ships for each young group we exclude known brown dwarfs. 
The reason for this is that the low-mass cut-off of the semi- 
empirical model isochrones we use occurs at ~ 0.1 Mq (see 
Section 4.1). Table [I] details our list of members for each 
group, along with the compiled VJ photometry, spectral 
types and distances. 















































Table 1. VJ photometry, spectral types and distances for member stars in our sample of young moving groups. The final column 
indicates whether the given star has a position in CMD spac e wh ich is commensurate with group membership based on the best-fitting 
model as calculated using the r 2 fitting statistic (see Section [~4. 2 | for details). The full table is available as Supporting Information with 
the online version of the paper and includes all members of the AB Dor moving group, Argus, BPMG, Carina, Columba, r] Cha, TWA, 
Tuc-Hor and 32 Ori. 
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2.1 Members 


Our list of member stars was assembled from literature mem- 


bership lists and includes the following studies: |Mamajek, 

Lawson, & Feigelson (1999), Luhman & Steeghs (2004), Lyo 

et al. (2004), Zuckerman & Song (2004), 

Mamajek (20051, 

Lopez-Santiago et al. (2006), Torres et al. 

12006), Mamajek 

(2007 

1, Lepine & Simon (2009), Schlieder, Lepine, & Simon 

(2010 

1, Kiss et al. (2011), Rodriguez et al. (2011 

, Zuck- 

erman et al. (2011), Schlieder, Lepine, & Simon 

2012al, 

Schlieder, Lepine, & Simon (2012b), Shkolnik et al. 

(2012), 

Binks & Jeffries (2014 

), Ducourant et al. (2014), Kraus et al. 

(|2014|), Riedel et al. l 

2014), and Shvonski et al. (in prepa- 


ration). In addition we include the ‘bona fide’ and high- 
probability (> 90 per cent) candidate members as defined 


(2014b). 


by Malo et al. (20131, Malo et al. (2014a I and Malo et al. 


The series of papers by jMalo et al.| have identified sev¬ 
eral hundred high-probability candidate members for sev¬ 
eral of the groups included in this study, however many of 
these still require additional measurements (e.g. Li equiva¬ 
lent width, radial velocity, trigonometric parallax measure¬ 
ment) to unambiguously assign final membership to a given 
group. Given that our aim here is to derive the ‘best’ rep¬ 
resentative isochronal ages for young moving groups in the 
solar neighbourhood, it is therefore critical that we minimise 
the number of contaminating interlopers i.e. candidate mem¬ 
bers which still have a questionable status. Therefore, for the 
inclusion of |Malo et al. high-probability candidate members, 
we require that such stars must have a measured radial ve¬ 
locity and/or trigonometric parallax which are/is consistent 
with membership for a given group. 

Not accounting for unresolved multiples, our sample in¬ 
cludes a total of 89 members of the AB Dor moving group, 27 
members of Argus, 97 members of the BPMG, 12 members 
of Carina, 50 members of Columba, 18 members of r/ Cha, 
189 members of Tuc-Hor, 30 members of TWA, and 14 mem¬ 
bers of 32 Ori. 


2.2 Photometry 


The stars in our list are spread over a large area on the sky, 
and although there have been dedicated searches for low- 
mass members in some of the groups studied here (see e.g. 
Torres et al.||20d6 1, there is a significant dearth of homoge¬ 
neous optical photometric coverage of these moving groups 
(especially in the M dwarf regime). This is not the case 
in the near-infrared (near-IR), however, where all member 
stars have counterparts in the Two-Micron All-Sky Survey 
(2MASS; Skrutskie et al.j2006) Point Source Catalog (PSC; 
Cutri et al.|2003). 

Although ages can be derived using near-IR CMDs, the 
loci become vertical with colours J — K s ~ 0.9 mag for 
2500 < T/ff 4000 K. So although stars become less lumi¬ 
nous as a function of time, the sequences for different groups 
are almost degenerate with age. Furthermore, at very young 
ages (< lOMyr), observations are further complicated by 
the presence of circumstellar material. Therefore, to derive 
ages for our sample of groups we supplement the near-IR 
photometry with F-band data. Note that due to increased 
contamination in the K a -band photometry as a result of cir¬ 


cumstellar material (especially for r/ Cha and TWA), we 
derive ages using the My, V — J CMD. 

Whilst many of the stars in our list of members have 
Hipparcos entries, this catalogue, unfortunately, does not 
quote explicit uncertainties on the F-band magnitude. We 
have therefore assembled F-band photometry from the 


following sources: Vyssotsky (1956), Hauck & Mermilliod 
(19981, Tycho-2 (Hpg et al. 2000) Fr transformed into John- 

son F using the relation of 

Mamajek, Meyer, & Liebert| 

(2006), Lawson et al. (2001 

1, Lawson et al. (2002), the 

All-Sky Automated Survey (ASAS; 

Pojmanski 2002), Lyo 

|et al. (|2004), Zuckerman & Song i 

2004), the Naval Ob- 


servatory Merged Astrometric Dataset (NOMAD; Zacharias 
et al. 2005|), Barrado y Navascues ( 2006|), Costa et al. (20061, 


Mcrmilliod (2006|, the Search for Associations Containing 


Young stars (SACY) sample (Torres et al. 20061, the All- 
Sky Compiled Catalogue of 2.5 million stars (ASCC-2.5; 
|Kharchenko fe Roeser| |2009| ), |Chauvin et al.| (|2010|), the 


Southern Proper Motion Catalog 4 (SPM4; Girard et al. 


20111, Lepine & Gaidos (20111, the USNO CCD Astrograph 
Catalog 4 (UCAC4; Zacharias et al.|2013 l, and (Riedel et al. 

( Mi] ). 

Whilst the majority of these studies provide associated 
uncertainties on the F-band magnitude, stars for which we 
have adopted NOMAD and SPM4 photometry, as well as 
photometry taken from the studies ofiZuckcrm an fe Song| 
(2004|, Barrado y Navascues (20061, Chauvin et al. (20101, 
and Lepine & Gaidos (20111 do not. In such cases we adopt 
a conservative F-band uncertainty of ay = 0.1 mag. Follow¬ 
ing Pecaut & Mamajek (20131, for SACY objects brighter 
than F = 12 mag we adopt an uncertainty of ay = 0.01 mag, 
whereas for fainter sources we again adopt the rather con¬ 
servative ay =0.1 mag. 

There are two issues which we must consider before 
placing a given star in the CMD, namely unresolved binary 
stars and low-quality 2MASS PSC photometry. Addressing 
the former, our list of members naturally includes a signif¬ 
icant number of binaries, and although for several of these 
both components are resolved in the F- and J-bands, for 
others, the optical and/or near-IR photometry is unresolved. 
If the two components are resolved in both the F- and J- 
bands, then we plot each component separately in the CMD. 
However, given that we have heterogeneous information con¬ 
cerning the binary population of all the co-moving groups 
(e.g. the lack of information concerning the swathe of new 
Tuc-Hor members from K raus et al.|2014| coupled with the 
fact that the model isochrones we use to derive ages include 
the effects of binary stars (see Section 4.21, then if either 
the F- and/or J-band photometry is unresolved we instead 
opt to plot the combined system measurement in the CMD. 
Note that the model isochrones do not account for higher or¬ 
der multiples (e.g. triples or quadruples), and so if we have 
information concerning the number of components for such 
systems (e.g. TWA 4 is a quadruple system) we then attempt 
to derive individual component colours and magnitudes fol¬ 
lowing the technique of |MermTl liod et al. |1992) based on 
the system colour and magnitude as well as the magnitude 
difference between the components in a given photometric 
bandpass. This process is explained fully in Appendix |A| and 
specifies the systems for which this technique was adopted. 
In the case of higher order multiple systems for which we 
do not have the required information to derive individual 







































































































































































































































Solar neighbourhood age scale 5 


component colours and magnitudes (e.g. the magnitude dif¬ 
ference between the various components) we simply plot the 
combined system measurement in the CMD. 

Secondly, to ensure we are using the best possible pho¬ 
tometry from the 2MASS PSC, we only use those objects for 
which the associated Qflg is ‘A’. For many of the brighter, 
early-type stars, however, the associated Qflg is ‘C’, ‘D’ or 
‘E’ i.e. either a magnitude was extracted but the associated 
uncertainty is prohibitively large or there were serious is¬ 
sues extracting a magnitude at all. For several cases, whilst 
the A' s -band photometry has an associated Qflg of ‘A’, the 
J-band photometry does not, and so we therefore use the 


main-sequence relation of Pecaut & Mamajek (2013) to in¬ 


fer a V— J colour for the star based on its V — K s colour. For 
a smaller number of cases all three 2MASS PSC bands have 
associated Qflgs which are not ‘A’, and in such instances we 
instead use the associated B — V colour (typically from|Mer- 


milliod 2006) to infer the V—J colour again using the Pecaut 


& Mamajek 


12013) relation. The stars which are affected 


as such are expected to have settled on the main-sequence 
(based on their spectral type) and so we can be confident 
that such a technique will yield colours which can be used 
to reliably place stars in the CMD. A full list of the stars 
for which we have adopted this boot-strapping of colours, 
in addition to objects with questionable memberships and 
those we have excluded from our age analysis, can be found 
in Appendix [Bl 


2.3 Distances 

We preferentially adopt trigonometric parallax measure¬ 
ments for assigning distances to each star, however in cases 
where these are not available or the associated uncertainties 
are prohibitively large, we adopt either a kinematic distance 
from the literature or derive one using the so-called ‘con¬ 
vergent point’ method (de Bruijne 1999 Mamajek 20051. In 
certain cases, the trigonometric parallaxes are possibly er¬ 
roneous (e.g. TWA 9A and TWA 9B where the Hipparcos 
parallax may be in error at the ~ 3 a level; see Pecaut & Ma¬ 
majek 20131, and therefore we adopt a kinematic distance 


instead. Of the young groups in our sample we require kine¬ 
matic distances for several members of the AB Dor moving 
group, BPMG, Tuc-Hor and TWA, for which we adopt the 


convergent point solutions of Barenfeld et al. (20131, Mama¬ 


jek & Bell (20141, Kraus et al. (20141 and Weinberger et al. 
(20131 respectively. 


For the two most distant groups in our sample - r\ 
Cha and 32 Ori - there are only a small number of stars 
with trigonometric parallax measurements (RECX 2 [p Cha] 
and RECX 8 [RS Cha] in p Cha and 32 Ori, HR 1807 and 
HD 35714 in 32 Ori). As opposed to deriving individual kine¬ 
matic distances for the other members in these groups, we 
instead adopt distances of 94.27± 1.18pc and 91.86T2.42pc 
for stars in p Cha and 32 Ori respectively based on the 
weighted average of stars within these clusters with trigono¬ 
metric parallaxes from the revised Hipparcos reduction of 


3 THE MODELS 

In this study we adopt four sets of semi-empirical model 
isochrones based on different interior models, namely the 
Dotter et al. (20081, Tognelli et al. ( 2011[ |, [Bressan et al.| 
(20121 and Baraffe et al. (20151 models (hereafter referred to 
as Dartmouth, Pisa, PARSEC and BHAC15 respectively^ 
Note that the Pisa models represent a customised set of 
interior models which span a much greater mass and age 


range than those available via the Pisa webpages (see Bell 
et al.| 2014 for details). In addition, we note that the semi- 
empirical PARSEC models are not the recent vl.2S models 


by Chen et al. (2014). Instead, they are based on the vl.l 
interior models (computed assuming a solar composition of 
Z = 0.0152) and were created using an identical process to 
the other semi-empirical model isochrones; a process briefly 
summarised below. 

It has been demonstrated that using atmospheric mod¬ 
els to transform theoretical interior models into CMD space 
results in model isochrones that do not match the shape 


of the observed locus, especially at low-masses (see e.g. Bell 
et al. 2012 Lodieu 20131. Even with the use of empirical BC 


T e g relations, the aforementioned discrepancy between the 
models and the data still exists, especially in optical band- 
passes (see e.g. Bell et al.|2014 1 . Hence, in Bell et al. (2013 


20141 we described a method of calculating additional em¬ 


pirical corrections to the theoretical BC-Teft relations pre¬ 


dicted by the BT-Settl atmospheric models (Allard et al. 


2011) so as to match the observed colours of low-mass stars 


in the Pleiades assuming an age of 130 Myr, distance mod¬ 
ulus of 5.63 mag and reddening of E(B — V) = 0.04 mag. 
These corrections were then combined with the theoretical 
dependence of the BCs on surface gravity from the atmo¬ 
spheric models to derive semi-empirical BC-Teff relations. 

In [Baraffe et al. ||2015| the authors suggest that the new 
BHAC15 models now match the shape of the observed lo¬ 
cus, in the sense that the optical colours for M dwarfs are no 
longer too blue. In Fig. [T]we show the V, V — J CMD of the 
Pleiades and overlay 130 Myr model isochrones from differ¬ 
ent ‘generations’ of the Baraffe et al. models (assuming the 
same distance modulus and reddening as noted in the previ¬ 
ous paragraph). From Fig. [l] it is clear that whilst progress 
has been made since the original Baraffe et al. (1998) models 


were first published, there are still unresolved problems with 
the latest BHAC15 in the optical and/or near-IR bandpasses 
at low Teff ■ Repeating the binary analysis of Bell et al. (2012) 


van Leeuwen (2007). 


in which we compared the theoretical system magnitudes 
predicted by the theoretical models with the measured sys¬ 
tem magnitudes in the VI C JHK S bandpasses, we find that, 
for the BHAC15 models, the observed spread in the optical 
bandpasses is smaller than for the other models, however it is 
still appreciably larger than in the near-IR bandpasses. Thus 
despite the fact that the discrepancy between observed loci 
and the models is smaller when using the BHAC15 models 
(compared with previous ‘generations’), it is still necessary 
to calculate additional empirical corrections to the theoret¬ 
ical BC-Tefr relations before using these to derive absolute 
ages for young stellar populations using CMDs. 


1 The semi-empirical Dartmouth and Pisa models we use here 
are available via the Cluster Collaboration isochrone server http: 
//www.astro.ex.ac.uk/people/timn/isochrones/ 
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Figure 1. V, V — J CMD of the Pleiades. Overlaid are different 
‘generations’ of the Baraffe et al. model isochrones with ages of 
130 Myr, each of which has been reddened by the equivalent of 
E(B — V ) = 0.04 mag and shifted vertically by a distance modulus 
dm = 5.63 mag. The original Baraffe et al. l | 1998^ models (with a 
solar-calibrated mixing-length parameter a = 1.9) coupled with 
the NextGen atmospheric models are shown as the green line. The 
same interior models, but coupled with the BT-Settl atmospheric 
models (computed using the |AspIun d et al.||200 9 abundances) 
are represented by the blue line. The latest Baraffe ct al. ( |2015[ ) 
interior models coupled with the BT-Settl models (both computed 
using the[ Caffau et al.|2011| abundances) are denoted by the red 
line, ft is clear that whilst the discrepancy between the models 
and the data has decreased with time, there are still unresolved 
problems in the optical and/or near-IR bandpasses at low T e g. 


4 ISOCHRONAL AGES FOR YOUNG, 
NEARBY GROUPS 


We use the maximum-likelihood r 2 fitting statistic of Naylor 


& Jeffries (20061 and Naylor (20091 to derive ages from the 


My, V — J CMDs of our sample of young moving groups. 
The t 2 fitting statistic can be viewed as a generalisation of 
the \ 2 statistic to two dimensions in which both the model 
isochrone and the uncertainties for each datapoint are two- 
dimensional distributions. Not only does the use of the r 2 
fitting statistic remove the issue of objectivity introduced 
through ‘by-eye’ fitting of model isochrones, but it also al¬ 
lows us to include the effects of binarity in our model dis¬ 
tribution, provides reliable uncertainties on the derived pa¬ 
rameters and allows us to test whether the model provides 
a good fit to the data. Similarly to \ 2 j the best-fit model is 
found by minimising r 2 . 


4.1 Creating the model distributions 

The two-dimensional model distributions are generated us¬ 
ing a Monte Carlo method to simulate 10 6 stars over a given 
mass range. To populate our models we adopt the canonical 
broken power law mass function of|Dabringhausen, Hilker, fc| 


Kroupa (2008). An important feature of our fitting method 


is that we include the effects of binary stars, thereby modify¬ 
ing the model isochrone at a given age from a curve in CMD 
space to a two-dimensional probability distribution. To in¬ 
clude binary stars in our fits we assume a binary fraction 
of 50 per cent and a uniform secondary distribution ranging 
from zero to the mass of the primary. 

Our sample of member stars ranges from early B-type 
stars to late M-type dwarfs, and given that across this range 
there is an obvious trend for the observed binary fraction to 
decrease as a function of decreasing spectral type (or mass; 
see e.g. Duchene fc Kraus|2013 l, our choice of a uniform bi¬ 
nary fraction of 50 per cent for all spectral types is somewhat 
idealised. The effect of varying the adopted binary fraction 
in model distributions has previously been investigated by 


Naylor & Jeffries (2006) who demonstrated that even at un¬ 


realistic values (e.g. 80 per cent across all spectral types) the 
best-fit age is affected at the < 10 per cent level. The model 
grids we calculate are dense and cover a large age range 
[log(age) = 6.0 — 10.0 in steps of Alog(age) = 0.01 dex], and 
so given this insensitivity to the adopted binary fraction, 
we adopt a uniform fraction of 50 per cent for the ease of 
computing them. 

We note that due to the processes involved in creat¬ 
ing the semi-empirical BC-T e ff relations (see Section [3}, the 
lower T e ff limit of the BCs may be reached before the lower 
mass limit of the interior models. Effectively this means that 
the semi-empirical pre-MS model isochrones only extend 
down to masses of ~ 0.1 Mq. As a result, simulated sec¬ 
ondary stars with masses corresponding to T e s values below 
this T/ff threshold are assumed to make a negligible contri¬ 
bution to the system flux i.e. this is equivalent to placing the 
binary on the single-star sequence. Essentially, this means 
that we begin to lose the binary population before the single¬ 
star population in our two-dimensional model distributions 
and results in a ‘binary wedge’ of zero probability between 
the binary and single-star sequences at low masses (see e.g. 
Fig# 


4.2 Fitting the CMDs 


Whilst we have made every effort to include only high- 
probability (> 90 per cent) candidate members in our sam¬ 
ple of young moving groups, it is likely that a fraction of 
these stars are in fact non-members. We have furthe r im- 
proved the r 2 fitting statistic since the modifications of |Nay-| 
lor (2009 1 and the reader is referred to Appendix [C in which 


we introduce an updated method (cf. Bell et al. 


2013} to 


deal with non-member contamination by assuming a uni¬ 
form non-member distribution. We choose such a concep¬ 
tual framework, in part because it is better suited to the 
problem in hand, but also because it allows us to calculate 
a goodness-of-fit parameter, a step which was missing from 


our earlier soft-clipping technique used in Bell et al. (20131. 


Prior to deriving ages from the My, V — J CMDs for our 
sample of young groups, we must first address whether we 
need to account for local levels of interstellar extinction and 
reddening in the solar neighbourhood. All of our member 
stars he within the so-called Local Bubble (< 100pc). A re¬ 
cent study by Reis et al. (2011) performed an analysis of the 


interstellar reddening in the Local Bubble using Stromgren 
photometry and provided estimates of how E(b—y) varies as 
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Figure 2. Best-fitting My,V — J CMDs of the AB Dor moving group. The red circles represent fitted data, whereas the blue squares 
denote objects which are removed prior to fitting (see Section |4.2| for details). Top left: BHAC15. Top right: Dartmouth. Bottom left: 
PARSEC. Bottom right: Pisa. The dashed line in each panel represents the ZAMS relation for that specific set of model isochrones. 


a function of distance. To calculate how this local reddening 
could affect the colours and magnitudes of our member stars, 
we assume that Av = 4.3 x E(b — y) and Aj/Ay = 0.282 
(see Glaspey|1971 and Rieke & Lcbofs ky|1985 respectively) 
to infer that the typical levels of extinction and reddening 
are Av — 0.03 mag and E(V — J) — 0.02 mag respectively. 
Compared to the combined photometric and distance uncer¬ 
tainties, such levels of extinction and reddening are insignif¬ 
icant and therefore we can leave the semi-empirical models 
in the absolute magnitude-intrinsic colour plane to derive 
ages from the My, V — J CMDs. 

To illustrate our fitting procedure we will use the AB 
Dor moving group as an example, for which the best-fitting 
My, V — J CMDs are shown in Fig. [2] Prior membership 
probabilities for the individual stars are taken from the anal¬ 
yses of |Malo et al.| ( [2013||2014a|b| >, and for this dataset the 
mean prior probability of a given star being a member of the 


cluster P(M C ) is greater than 0.99 (see Appendix Cl.2 for an 
explanation of how the r 2 fitting statistic can be used to in¬ 
clude prior membership probabilities for individual sources). 
Therefore, with 89 stars in our catalogue, this would then 
imply that, at most, one object does not originate from the 
model sequence. Note that for member stars in our sample 
of young moving groups which are not included in the anal¬ 
yses of|Malo et al.|we adopt a uniform prior probability of 
0.9. 

An examination of Fig. [2j however, suggests that such 
a low non-member fraction is unlikely. There is a group of 
four objects which lie significantly above the equal-mass 


binary envelope of the model distribution between 1.5 < 
V — J < 2.5 mag, the brightest of which is brighter than 
an equal-mass quintuple. Furthermore, there are several ad¬ 
ditional objects which, even allowing for the uncertainties 
in My, also lie above the upper envelope of the best-fit 
model. It is worth noting that none of these objects are sin¬ 
gle stars, but are either resolved components of binary sys¬ 
tems or unresolved binaries, all of which have separations of 
< 4 arcsec, and therefore the effects of orbital motion could 
conceivably affect the parallax solutions of |van Leeuwen| 
12007). Note that there is also a faint object (GJ 393) at 


My ~ 10.4 mag, V — J — 3.5 mag, which is a single star 
with a well-constrained distance, that lies ~ 0.5 mag below 
the single-star sequence. Such objects therefore imply that 
the probability a given star is not modelled by our sequence 
is greater than the 0.01 probability suggested by the raw 
values of P(M C ). This conclusion is justified by examining 
the distribution of the r 2 values for each individual data- 
point. The top left panel of Fig. [3] shows the cumulative 
distribution for these values, along with the predicted dis¬ 
tribution calculated from the same 1000 simulated clusters 
(in this case based on the Dartmouth model distributions) 
used to calculate the goodness-of-fit parameter [Pr(r 2 ); see 
Appendix C3.1| . There is an obvious tail of points with val¬ 
ues much higher than the prediction from the simulations. 

To overcome this problem we therefore set a maximum 
prior membership probability (P max ) and multiply all the 
prior membership probabilities by this factor. We assessed 
the correct value for P max by adjusting it until the total of 
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Figure 3. Comparing the real and expected r 2 distributions (us¬ 
ing the Dartmouth models) for the AB Dor moving group. The 
black histograms show the fraction of stars with individual r 2 
values less than a given value in the fit to the data. The red 
curves represent the expected distributions of the probability of 
obtaining a given t 2 calculated for the simulated observation (us¬ 
ing the best-fit model). Each panel is for a different value of the 
maximum allowed prior membership probability P max (see text). 


the prior and posterior memberships roughly matched. This 
process is illustrated in Fig. [3] The drop at the end of the 
distribution is the stars which are apparent non-members, 
and this is fairly closely matched at P max = 0.94 which is 
the value where the total prior and posterior memberships 
are best matched. For the AB Dor moving group we find 
that, regardless of which set of semi-empirical pre-MS mod¬ 
els is adopted, the maximum prior membership probability 
is P max > 0.9, and that this results in goodness-of-fit values 
of Pr(r 2 ) > 0.5 in all cases. 

An important feature to bear in mind when calculat¬ 
ing reliable best-fit ages and goodness-of-fit parameters is 
that there are certain regions in CMD space which our grid 
of models [covering the age range log(age) = 6.0 — 10 . 0 ] 
simply do not occupy. Note that this is compounded for 
the BHAC15 models which have an upper mass cut-off of 
1.4 Mq, which means that any ‘handle’ on the age provided 
by higher mass stars (specifically those evolving between the 
zero-age main-sequence [ZAMS] and the terminal-age main- 
sequence [TAMS]) is effectively ignored by these models, 
which is not the case for the other models which have upper 
mass cut-offs of > 5 Mq . Although we have tried to exclude 
known brown dwarfs from our membership lists, there are 
still some stars which occupy regions of the CMD not cov¬ 
ered by the models (at both high- and low-masses; see e.g. 
the blue squares in the upper right panel of Fig. [9j, which 
must be removed from the fit before calculating either of the 
aforementioned parameters. 

A further consequence of our fitting technique is that 
our best-fit model returns posterior membership probabili¬ 
ties. For a given group in our sample we expect a negligible 
age spread (or equivalently luminosity spread), and there¬ 
fore we can use the posteriors to identify stars which appear 
to be non-members (based solely on CMD position) in an 
effort to further refine the membership lists of these young 
moving groups. Such a technique is far from unassailable. 


For example, we need only think of the obscuring effects 
of edge-on discs or uncertainties in the distances to objects 
(see e.g. TWA 9 in Section 2.31 which will act to shift stars 
away from the observed locus of the group. Identification of 
such stars, however, represents a sensible ‘first-order sanity 
check’ which should be used in combination with other met¬ 
rics (e.g. Li equivalent widths, radial velocities, etc.) to es¬ 
tablish whether a given star is a high-probability candidate 
member of a given moving group. For example, the anal¬ 
yses of |Malo et al.| have identified several high-probability 
candidate members in the BPMG, however an examination 
of Fig. [5] clearly shows that the inclusion of these stars re¬ 
sults in a luminosity spread of ~ 3 mag in My at a colour of 
V — J — 4.5 mag. Such a spread, given its age of > 20 Myr, is 
incomprehensibly vast and the fact that the faintest stars in 
this sample lie below the computed ZAMS lends confidence 
to our assertion that our membership list still contains likely 
non-members. This is also borne out by the calculated pos¬ 
terior probabilities for these stars which are all less than 0.1 
i.e. indicating that these are likely non-members. 

Figs. [4] |TTj show the best-fitting My,V — J CMDs for 
the remainder of our sample of young groups. Of the re¬ 
maining groups, we find that, except for Argus, the BPMG 
and 77 Cha, the maximum prior membership probability is 
P max > 0.9 and that the resultant goodness-of-fit values 
Pr(r 2 ) are >0.5 for all four sets of semi-empirical pre-MS 
models. It is likely that the lower P max values required for 
the BPMG and 77 Cha (of between 0.7 and 0.8) are a re¬ 
sult of contamination by non-members in the former (see 
above) and astrophysical phenomena affecting the CMD po¬ 
sitions of ~ 10 per cent of the members in the latter (see 
Section 5.51 in conjunction with our adoption of a well- 
constrained distance (< 2 per cent uncertainty) for all mem¬ 
ber stars. Whilst the derived P max values for Argus are all 
> 0.85, the corresponding goodness-of-fit values Pr(r 2 ) are 
< 0.2 for all of the adopted models. Again, it is likely that 
this is suggestive of a ~ 10—15 per cent level contamination 
in our membership list for this association (see Section 5.31. 

Table [2] shows the best-fit age for each group according 
to the four sets of semi-empirical models, in addition to our 
final adopted age. For the final age we adopt the average of 


i) the median, ii) the Chauvenet clipped mean (Bevington 
fc Robinson|1992|), and iii) the probit mean (|Lutz fe Upgren 


1980) of the individual best-fit ages. The quoted uncertain¬ 


ties on the final ages represent the statistical and system¬ 
atic uncertainties added in quadrature. Given the asymmet¬ 
ric statistical uncertainties on several of the individual ages 
shown in Table [2] we calculate both an upper and lower sta¬ 
tistical uncertainty by taking the median of the four individ¬ 
ual upper and lower statistical uncertainties. Our estimate 
of the systematic uncertainty arising from the use of differ¬ 
ent model isochrones is based on the average of i) the 68 per 
cent confidence levels and ii) the probit standard deviation 
I Lutz & Upgren]T980) of the individual best-fit ages. 


5 DISCUSSION 

In Section |4.2| we used semi-empirical models to derive a 
self-consistent, absolute isochronal age scale for nine young, 
nearby moving groups within 100 pc of the Sun. To assess 
the reliability of these ages, we must first compare these to 
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Figure 4. Best-fitting My,V — J CMDs of Argus. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 




V-J [mag] V-J [mag] 



Figure 5. Best-fitting My,V — J CMDs of the BPMG. The coloured symbols and dashed lines are the same as those in Fig. [2] Top 
left: BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 
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Figure 6. Best-fitting My, V — J CMDs of Carina. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 




Figure 7. Best-fitting My, V — J CMDs of Columba. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 
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Figure 8. Best-fitting My ,V — J CMDs of rj Cha. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 
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Figure 9. Best-fitting My, V — J CMDs of Tuc-Hor. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 
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V-J [mag] V-J [mag] 

Figure 10. Best-fitting My , V — J CMDs of TWA. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 




Figure 11. Best-fitting My, V — J CMDs of 32 Ori. The coloured symbols and dashed lines are the same as those in Fig. [2] Top left: 
BHAC15. Top right: Dartmouth. Bottom left: PARSEC. Bottom right: Pisa. 
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Table 2. Ages for the young groups studied in this paper. Individual ages are shown for each set of semi-empirical pre-MS model 
isochroness and have been derived using the r 2 fitting statistic for which we have set P max so that the total prior and posterior 
memberships are roughly equal (see Section |4.2| for details). The penultimate row lists our final adopted age for each group on which 
the quoted uncertainties represent the statistical and systematic uncertainties added in quadrature. The final row shows literature LDB 
ages for the BPMG and Tuc-Hor (see Section|5.1|for references) which highlights the consistency between the two age diagnostics. 


Model 




Group age (Myr) 






AB Dor 

Argus a 

BPMG 

Carina 

Columba 

ri Cha 

Tuc-Hor 

TWA 

32 Ori 

BHAC15 

145+889 

69^^® 

25 ± 1 

46± 46 

44±® 

12 ± 1 

46±* 

10 ± 1 

23t|| 

Dartmouth 

135 + g 5 

60l® 4 

23 ±1 

49l 43 

43t 8 3 

10 ± 1 

48± 4 

7+2 
'-1 

20±® 

PARSEC 

151 tli 

60+268 


45tL 

43± 3 

14 ± 1 

46 ± 1 

13 ± 1 

25ls 

Pisa 

166±m 

55tf 5 

20t^ 

00 50 

+ 1 

38 ±3 

8 ± 1 

40± 4 

9 ± 1 

18 ± 1 

Adopted 

149t“ 

- 

24±3 

45± 41 

CO 

+ 1 

CN 

11±3 

45±4 

10±3 

22t 4 

LDB age 

- 

- 

24 ±5 

- 

- 

- 

40 ± 3 

- 

- 


Notes: 

“We do not provide a final adopted age for Argus as it remains unclear whether the stars in our list of members are representative of a 
single population of coeval stars (see Section |5.3| for details). Note also that the best-fit age for Argus using the PARSEC models only 
provides an upper limit on the age of the association and hence we do not include a lower age uncertainty in the table. 


what are considered to be well-constrained ages for the same 
groups. 


ages against. Hence, we must place the isochronal ages for 
the seven remaining groups in context by comparing them 
to other age-dating techniques adopted in the literature. 


5.1 Comparison to LDB ages 

As discussed in the Introduction, LDB ages are arguably the 
most reliable age diagnostic - in terms of calculating abso¬ 
lute ages - currently available for stellar populations with 
ages of between ~ 20 and 200 Myr. Note that although LDB 
ages have been advocated as our best hope of establishing a 
reliable age scale in this age range (primarily due to the fact 
that LDB ages have been shown to be relatively insensitive 
to variations in the physical inputs adopted in the evolution¬ 
ary models and the fact that between different models there 


is excellent agreement; see e.g. Soderblom et al. 20141, recent 


studies have demonstrated that by not accounting for the ef¬ 
fects of starspots in the evolution of pre-MS stars, the LDB 
age scale may in fact only be good to ~ 10 — 20 per cent (see 


e.g. Jackson fc Jeffries|2014| |Somers fe Pinsonneault 2015[). 


Binks & Jeffries ( |2014[ ) and Malo et al. (2014b I have 
both identified the LDB in the BPMG, deriving ages of 21±4 
and 26 ± 3 Myr respectively (the difference in age stemming 
from the use of stellar evolutionary models which include 
a magnetic field model in the latter study). Despite this 
difference, the derived LDB ages are consistent and for the 
purposes of comparing this to our isochronal age we adopt 
a value of 24 ± 5 Myr. In addition to the available LDB age 


for the BPMG, Kraus et al. (2014) recently identified the 


LDB in Tuc-Hor, for which we adopt an age of 40 ± 3 Myr. 
It is therefore reassuring then to find that the isochronal 
ages we have calculated for both the BPMG and Tuc-Hor 
are consistent with the derived LDB ages. Note that our age 
for the BPMG also compares favourably to the 22 ± 3 Myr 


derived by Mamajek & Bell (20141 which was based on an 


My, B — V CMD analysis of the A-, F- and G-type members 
of the group. Whilst such consistency lends confidence to 
our absolute age scale for the young moving groups in the 
solar neighbourhood, there are (currently) no other LDB 
ages for our sample with which to compare the isochronal 


5.2 AB Dor moving group 

Literature ages for the AB Dor moving group vary from 


the relatively young (~ 50 — 70 Myr; see e.g. Zuckerman, 


Song, & Bessell 2004; To rres et al.|2008 |da Silva et al.|2009 [ 
to essentially coeval with the Pleiades (~ 130 Myr; see e.g. 
Luhman, Stauffer, & Mamajek 2005). A more recent analysis 


by Barenfeld et al. (20131 provided a strong constraint on 


the age of the group by identifying the main-sequence turn¬ 
on for AB Dor nucleus stars in the My, V — K s CMD. They 
demonstrated that the late K-type stars have already settled 
onto the ZAMS, and used this to place a firm lower limit on 
the age of 110 Myr. 

Interestingly, the recent discovery of a Li-rich M8 
high-probability brown dwarf candidate member (2MASS 
J00192626+4614078; |Gagne et al. 2014 1 provides us with 
a strong upper limit on the age of the group of 196 Myr 
(Binks, private communication). Additional candidate mem¬ 
bers with spectral types in the range M4-M8 have been re¬ 


ported by Gagne et al. (2015) and thus spectroscopic follow¬ 


up of these objects could place stronger constraints on the 
absolute age of the AB Dor moving group. Our isochronal 
age therefore lies almost exactly midway between the lower 
and upper age limits provided by the main-sequence turn¬ 
on and the lower edge of the tentative identification of the 
LDB. 

From Table [2j it is clear that the main ‘handle’ on the 
age of the AB Dor moving group comes from the high-mass 
stars which are evolving between the ZAMS and TAMS. 
The BHAC15 models demonstrate that if these stars are 
omitted from the fit then, although the best-fit age is in 
general agreement with those from the other models, the as¬ 
sociated upper uncertainty on the age is prohibitively large. 
The reason why the low-mass population of the group can¬ 
not provide a precise age is due simply to the fact that at 
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ages of > 100 Myr the model isochrones occupy essentially 
the same position in CMD space (except for very-low-mass 
objects with V — J > 4.5 mag for which we have only a 
few stars in our AB Dor moving group membership list; see 
Fig.§, resulting in a relatively flat r 2 space across a wide 
range of ages. 


5.3 Argus association 

The kinematics of Argus are very similar to those of the 
nearby young cluster IC 2391 (Torres et al.|2003|, for which 


Barrado y Navascues, Stauffer, & Jayawardhana (2004) de¬ 


rived an LDB age of 50 ± 5 Myr. A similar age for Argus 
(~ 40 Myr) was proposed by Torres et al. (20081 on the basis 


of Li equivalent widths and positions in the V, V — I c CMD, 
and therefore our isochronal ages of ~ 60 Myr disagree with 
the current literature age for this association. Furthermore, 
whilst the different sets of semi-empirical models tend to 
agree upon an age of ~ 60 Myr for Argus, the uncertain¬ 
ties on these individual ages are extremely large and imply 
that based on our current membership list, it is not valid to 
assign a unique, unambiguous age to Argus (see Table [2]). 
We believe that the smaller uncertainties on the BHAC15 
best-fit age stem from the upper-mass cut-off of these mod¬ 
els which neglect any age information from the higher mass 
members. Excluding the BHAC15 models for the moment, 
the other ages in Table [2] clearly demonstrate that although 
there is a minimum r 2 value within the grid, the r 2 space 
is also relatively flat over a wide range of ages (see e.g. the 
PARSEC models which can only give an upper age limit of 
328 Myr). 

Our membership list for Argus stems primarily from the 
Bayesian analyses of Malo et al. ( |2013 2014a|b I and com¬ 
paring it with the membership list of Torres et al.| ( |2008| 
their Table 12) there are two notable differences i) only three 
stars in common appear in both lists (BW Phe [HD 5578], 
HD 84075 and NY Aps [HD 133813]) and ii) the median 
distance of members in our list is 32 pc whereas for the |Tor-| 
|res et aT~| list it is 96 pc. Furthermore, looking also at the 
distribution of distances (from the nearest to the farthest) 
in these two membership lists, those stars in our list cover a 


range of ~ 65 pc, whereas those in Torres et al. (20081 cover 


a range of more than twice this at ~ 150 pc (cf. ~ 50 pc 
for the supposedly coeval Tuc-Hor). Given this spread in 
distances for the |Torres et al .] members, one can reasonably 
ask whether we expect these objects to have formed together 
and therefore be coeval in the first place. 

Fig.0 shows the best-fitting CMDs of Argus and there 
are two points worth mentioning. First, there is a group of 
5 A-type stars (all of which were proposed as members by 
IZuckerman et al.l2011l and none of which are unresolved bi¬ 
naries [according to the Washington Visual Double Star Cat¬ 
alog; jMason^rjal.][200l]), which we would expect to define 
the ZAMS of the association (as is the case with the other 
groups in our sample). None of these, however, actually lie 
on the ZAMS (4 are over-luminous and 1 is under-luminous 
with respect to the ZAMS), which suggests that these stars 
may not be coeval, but rather represent stars at different 
evolutionary stages. Second, of the K- and M-type stars in 
the association, a significant fraction of these appear to lie 
below the best-fitting model in all panels of Fig. [4] Whilst 
these stars all appear to be young and active, as evidenced 


by a combination of high Rx = Lx/Lbol values and large 
Ha equivalent widths (see e.g. Riaz, Gizis, & Harvin|20061 , 
they all lack Li measurements which may help us discrimi¬ 
nate between young active stars belonging to a given group 
and slightly older active ZAMS stars which do not (at least 
for the mid to late M-type objects). Given these ambigui¬ 
ties, we are therefore reluctant (at present) to assign a final, 
unambiguous age to Argus as it appears as though our mem¬ 
bership list suffers from a high level of contamination, and 
hence it remains unclear whether this list represents a single, 
coeval population of stars, or even whether the association 
is in fact physical. 


5.4 Carina and Columba associations 

Carina is an extremely sparse association whose Galactic 
space motion is, to within the uncertainties, statistically in¬ 


distinguishable from that of Columba (see e.g. Malo et al. 


2014aI, although spatially the two associations are rather 


distinct. Both associations were identified by |Torres et al.| 
12008) who further demonstrated that they share a simi¬ 
lar age of ~ 40 Myr through a combination of Li equivalent 
widths and positions in the V, V — I c CMD. Our isochronal 
ages for both Carina and Columba are consistent with ages 
of ~ 40 Myr and, to within the uncertainties, appear to be 
coeval (as well as share a common age with Tuc-Hor; see 
Torres et al. |2008 ). 

WE demonstrates that our best-fit models appear to 
follow the observed locus of Columba (tracing both the up¬ 
per and lower envelopes of the relatively well-populated as¬ 
sociation). For the much sparser Carina association, how¬ 
ever, such agreement between the models and the data 
is not so obvious [see Fig. [6] this is also reflected in the 
systematically lower Pr(r 2 ) values for Carina compared to 
Columba], Although we only have 12 stars in our mem¬ 
bership list for Carina, there is a case to be made that 
there are two apparent outliers which appear to lie below 
the observed locus. These stars (namely 2MASS J04082685- 
7844471 and 2MASS J09032434-6348330) will act to ‘drag’ 
the best-fit model to older ages, the result of which, is that 
2MASS J10140807-7636327 (the reddest star in our list) 
is assigned a posterior membership probability which indi¬ 
cates non-member status, despite occupying a position in 
CMD space which is commensurate with its binary status 
i.e. ~ 0.5 mag above the single-star sequence. 

If we repeat our isochronal age analysis having re¬ 
moved 2MASS J04082685-7844471 and 2MASS J09032434- 
6348330, we find a slightly younger age for Carina of ~ 
36 Myr (cf. 45 Myr). It remains unclear whether the kine¬ 
matic distance estimates for these two stars are erroneous 
or whether they are older than the other stars in our list. 
To better constrain the age of Carina, it is therefore clear 
that we have to determine whether these stars are genuine 
members of Carina or whether they are simply older ac¬ 
tive stars (akin to what we previously discussed in Argus; 


see Section 5.31. Given the spectral types of these two stars 


(early M-type) it is unclear whether further spectroscopic 
information (e.g. Li equivalent widths) would allow us to 
differentiate between these two options, and so this may be 
a case of having to wait until Gaia provides us with the 
necessary kinematic information to unambiguously do so. 
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5.5 r] Cha cluster 

The use of the LDB technique to derive ages is restricted to 
populations with ages of > 20Myr and thus ri Cha is too 
young for the adoption of such a technique. There are, how¬ 
ever, model-independent methods of deriving an age which 
can then be used as an independent diagnostic to compare 
against isochronal ages. One such method involves calculat¬ 
ing the time of minimum separation between stellar groups 
in the past on the assumption that the groups share a com¬ 
mon origin. As with other kinematic methods this approach 
has led to contradictory conclusions. For example, | Jilinski,] 


Ortega, & de la Reza (2005) performed a kinematic trace- 


back of the i] Cha and e Cha clusters assuming a Galactic 
potential and found that the smallest separation of only a 
few pc was ~ 7Myr ago. This coevality of the two clusters 


was later refuted by Lawson et al. (2009) who demonstrated 


that both clusters have significantly different ages based on 
a combination of CMD positions and surface gravity indi¬ 
cators. Furthermore, a more recent kinematic analysis by 
|Murphy~t al. (j2013) concluded that not only is r/ Cha a 
few Myr older than e Cha, but also that there is little ev¬ 
idence that the two clusters were appreciably closer than 
their current ~ 30 pc in the past. 

The literature isochronal ages for 77 Cha imply an age 
of 5 — 8 Myr (e.g. Luhman & Steeghs|[5004|, however this 
becomes slightly older when the effects of binary stars are 
accounted for (cf. 7—9 Myr; Lyo et al. 2004). The fundamen¬ 
tal parameters of the double-lined eclipsing binary RECX 8 
(RS Cha) are arguably some of the best constrained in our 
entire sample of member stars, with typical uncertainties 
of only a few per cent. Such precise measurements naturally 
provide stringent tests of stellar evolutionary models by forc¬ 
ing them to fit both components of the binary system at a 
given age (see e.g. |Stassun, Feiden, fe Torres|[20 14l. Simul¬ 
taneous fitting of the masses, radii, T e fr and Lboi for both 
components by |Alecian et al. (20071 and Gennaro, Prada 


|Moroni, fe Tognelli| ( |2012| ) imply an age of ~ 9 Myr for the 
RECX 8 system, which is consistent with our isochronal 
age of ~ 11 Myr for the cluster. Furthermore, Fig. [ 8 ] shows 
that our semi-empirical pre-MS models provide a good fit 
to the data, tracing both the lower and upper envelope of 
the observed locus, except for the two outliers RECX 13 and 
RECX 15. 


Luhman & Steeghs (2004) demonstrated that unlike the 


other members of Cha, which exhibit negligible levels of 
interstellar extinction, RECX 13 has a measured extinction 
of Av = 0.4 mag [which corresponds to E(V—J) = 0.29 mag 
based on the Aj/Av = 0.282 relation of|Rieke fe Lebofsky| 


1985 see Section |4.2|. De-reddening RECX 13 using these 


values would place the star much closer to the single-star 
sequence of the model distribution. RECX 15, on the other 
hand, has a large IR excess at wavelengths of > 2 //in which 


is indicative of circumstellar material (see e.g. Megeath et al. 
|2005[ ). The inclination of this disc with respect to our line- 
of-sight remains ill-constrained, however high inclinations 
for discs around other r] Cha members have been reported 
(> 60°; |Lawson, Lyo, & Muzerolle 2004|. If the disc around 
RECX 15 is so inclined, then we would expect significant 
dimming in the optical wavelengths as a result of observing 
the star through it’s disc, and this could therefore explain 
why RECX 15 appears fainter than the other stars which 


represent the lower envelope of the cluster locus in CMD 
space. 


5.6 TWA 

TWA was the first of the young moving groups in our sam¬ 


ple to be identified in the literature (see e.g. de la Reza et al. 
|1989[ |Kastner et al.|[l997[ ) and as such many age estimates 
are now available for this association. Several authors have 
attempted to use kinematic information to derive a model- 
independent age for the TWA, however, they are either con¬ 
tradictory or cover a prohibitively large range to provide a 


strong constraint on the age. For example, Mamajek (2005) 


noted that whilst the data are consistent with expansion, the 
corresponding expansion age of 20+y 5 Myr has such large as¬ 
sociated uncertainties that it is of limited use. 

More recent analyses by [Weinberger, Anglada-E scude(1 


& Boss (2013) and Ducourant et al. (20141, which have mea¬ 


sured a larger number of parallaxes and proper motions than 


in the Mamajek (2005) study, have come to completely op¬ 


posite conclusions. Weinberger et al. (20131 find that the 


space motions of TWA members are essentially parallel and 
do not indicate convergence at any time in the past 15 Myr, 


whereas Ducourant et al. (20141 find that a subset of 16 
members occupied the smallest volume of space ~ 7.5 Myr 


ago. Note that in the Ducourant et al. (2014) study, of the 


25 stars with reliable radial velocity and parallax informa¬ 
tion, the authors removed 9 stars from the sample (36 per 
cent) prior to deriving the traceback age as they demon¬ 
strated discrepant space motions. Omissions such as these 
again raise the question of subjectivity concerning the stars 
which one includes/excludes in a given sample. 

Despite the lack of an unambiguous kinematic age for 
the TWA, several isochronal ages have been reported in 
the literature which suggest an age of ~ 10 Myr (see e.g. 
Webb et al.||1999l |Barrado y Navascues||2006] |Weinberger| 


et aL|20T3l. Whilst our isochronal age is in excellent agree¬ 


ment with previous isochronal age estimates, Fig. |10| shows 
that the Dartmouth models (in particular) do not provide a 
good fit to the data. Whilst it is apparent that these mod¬ 
els have an enlarged ‘binary wedge’ when compared to the 
other models (presumably stemming from the lower-mass 
cut-off in the interior models and the slightly different T e fi 
scale) which could play a role in this discrepancy, it is also 
notable that there is a mismatch between the observed slope 
of the locus and that of the best-fitting model. Furthermore, 
even the other models appear to lie slightly above the lower 
envelope of the TWA locus. Hence, it is possible that our 
isochronal age could be underestimated and that a more 
representative age would be closer to ~ 15 Myr. 

Even when using models which include the effects of bi¬ 
nary stars, it is clear that we are unable to fit both the lower 
and upper envelopes of the observed locus in the CMD of 
the TWA. One reason for this increased luminosity spread 
could be due to obscuration effects arising from discs around 
the stars which define the lower envelope (namely TWA 6 , 
TWA 9B, TWA 21 and 2MASS J10252092-4241539) which, 
depending on their orientation, may act to make them ap¬ 


pear dimmer in the CMD (see e.g. RECX 15 in Section 5.51. 
A recent analysis of IR excess disc emission by Schneider, 
Melis, & Song (2012a I suggests that none of these stars have 
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appreciable excess emission and therefore this is unlikely to 
be the cause of the apparent luminosity spread. 

An alternative possibility is that the apparent luminos¬ 
ity spread could be a consequence of a real age spread within 
the association. Weinber ger et al7| ( |2013[ > propose that an age 
spread is probable given that the width of the age distribu¬ 
tion for TWA members exceeds that which would be ex¬ 
pected for a population of stars formed in a single burst. To 
place stricter constraints on whether the data are consistent 
with a real age spread or not, a more comprehensive study, 
such as that of .Reggiani et ah] (2011) in the Orion Nebula 
Cluster which accounted for uncertainties in the distance, 
spectral type, unresolved binarity, accretion and photomet¬ 
ric variability, would have to be performed. 


5.7 32 Ori group 


The 32 Ori group is fairly new and has not been as well- 
characterised as the other groups in this study. As such, 
our isochronal age represents the first definitive published 
age for the group. The group is located in northern Orion, 
centred near 32 Ori and Bcllatrix (although not containing 
the latter), with a distinctly high proper motion (p a , ps — 
+7, — 33masyr _1 ) and radial velocity ( v r — — 18.5kms _1 ) 
compared to the much more distant Ori OBI young stel¬ 
lar population. It was first reported by Mamajek (2007), 
and a Spitzer IR survey of the group reported by |Shvon-| 
[sk i et ah] ( 2010 ]) found several members to have dusty de¬ 
bris discs cooler than ~ 200 K, namely HD 35499 (F4V), 
HD 36338 (F4.5V), HR 1807 (AOVn), TYC 112-917-1 ( K4), 
and TYC 112-1486-1 (K3). |Kharchenko et al.| ( |2013[ ) re¬ 


covered the cluster and estimated it to to be at 95 pc, 
with mean proper motion p a , ps — +10.0, —32.2masyr -1 
(± 0 . 8 masyr _1 ), a core radius of ~ 1 . 6 pc, and an isochronal 
age of 32 Myr. 

Our age for 32 Ori is statistically indistinguishable from 
that for the BPMG (22 and 24 Myr respectively) and the 
coevality of these two groups is also evident from the po¬ 
sition of the main-sequence turn-on. In the 32 Ori group, 
HD 36338 (spectral type F4.5V) lies in a region of the CMD 
which pMamajek &; Bell| ( |2014[ ) termed the ‘false ZAMS’ i.e. 
in a region between the penultimate luminosity minima and 
the final luminosity minima corresponding to the ZAMS. All 
stars cooler than HD 36338 are over-luminous with respect 
to both the ‘false ZAMS’ and ZAMS, and therefore appear 
to be pre-MS stars. Note that the spectral type at which 
the main-sequence turn-on occurs is slightly later in 32 Ori 
than the BPMG (in which it occurs approximately one spec¬ 
tral type earlier) despite the slightly younger isochronal age 
for the group, however the ±3 Myr uncertainty on both 
isochronal ages and the paucity of stars with spectral types 
between F4.5 and F9 in 32 Ori makes such a direct compar¬ 
ison of only limited validity. A comprehensive characterisa¬ 
tion of both the stellar and circumstellar disc content of this 
group is currently underway (Shvonski et al. in preparation) 
which should help to better constrain this important epoch 
for terrestrial planet formation. 


6 CONCLUSIONS 

In this study we present a self-consistent, absolute age scale 
for eight young (< 200 Myr), nearby (< 100 pc) moving 
groups in the solar neighbourhood. We use previous liter¬ 
ature assessments to compile a list of member and high- 
probability candidate members for each of the young groups 
in our sample. Creating four sets of semi-empirical pre-MS 
model isochrones based on the observed colours of young 
stars in the Pleiades, in conjunction with theoretical cor¬ 
rections for the dependence on logg, we combine these with 
the t 2 maximum-likelihood fitting statistic to derive ages for 
each group using the My, V — J CMD. Our final adopted 
ages for each group in our sample are: 1491®], Myr for the AB 
Dor moving group, 24 + 3 Myr for the /3 Pic moving group 
(BPMG), 45ly J Myr for the Carina association, 42+| Myr 
for the Columba association, 11±3 Myr for the 77 Cha cluster, 
45 ± 4 Myr for the Tucana-Horologium moving group (Tuc- 
Hor), 10±3 Myr for the TW Hya association, and 22+g Myr 
for the 32 Ori group. At this stage, we are uncomfortable 
assigning a final, unambiguous age to Argus as it appears as 
though the stars in our membership list for the association 
suffer from a high level of contamination and hence may 
not represent a single, coeval population. Comparing our 
isochronal ages to literature LDB ages, which are currently 
only available for the BPMG and Tuc.-Hor, we find consis¬ 
tency between these two age diagnostics for both groups. 
This consistency instills conhdence that our self-consistent, 
absolute age scale for young, nearby moving groups is robust 
and hence we suggest that these ages be adopted for future 
studies of these groups. 
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APPENDIX A: BETWEEN A ROCK AND A 
HARD PLACE: DE-CONSTRUCTING 
PHOTOMETRY OF HIGHER-ORDER 
MULTIPLE SYSTEMS 


A1 Methodology 


To de-construct photometry of higher-order multiple sys¬ 
tems into individual components measurements we adopt 


the technique of Mermilliod et al. (1992 I which, in our case, 


uses the system V-band magnitude, system V — J colour 
and AV-band magnitude between the two components. 

We illustrate this technique using the TWA member 
TWA 4 as an example. TWA 4 is a quadruple system (un¬ 
resolved spectral type of K6IVe; Pecaut & Mamajek) 2013) 
comprising two visual components (separation < larcsec), 
which themselves are binary systems: A is a single-lined 
spectroscopic binary and B is a double-lined spectroscopic 
binary. Of the four components (designated Aa, Ab, Ba 
and Bb), only light from three of these - Aa, Ba and Bb 
- is visible in the combined optical spectrum (Soderblom 
et al. 1998). We adopt the system V-band magnitude of 


Vas = 8.912 ± 0.053 mag based on 601 observations by the 
ASAS. The 2MASS PSC gives a system J-band magnitude 
of Jab = 6.397 ±0.020 mag, yielding a system V — J colour 












20 C. P. M. Bell et al. 


of (V — J) ab = 2.515 ± 0.057mag. For the difference in 
V-band magnitude between the A and B components we 
adopt AVa.b = 0.504 ± 0.030 mag from |Soderblom et al.| 
(19981 based on two separate V-band measurements. 


Whilst the Mermilliod et al. (1992 I formalism is demon¬ 


strated for splitting BV photometry, this can be generalised 
to any combination of photometric bandpasses so long as 
the relationship between the colour and the magnitude is 
known (in our case the trend between My and V — J). We 
therefore fit an empirical (linear) slope to the resolved TWA 
members in the My, V — J CMD, for which we calculate a 
slope a = 1.530, and derive the colour and magnitude of the 
primary component according to 


(V-J) a = (V-J) AB + 2.5l0g 
and 


l _|_ jQ-0.4(4m)(l+«)/a 
^ _|_ JQ-0.4Am 


(Al) 


V A = Vab + 2.5 log(l + i 0 -°' 4Am ) (A2) 

where Am is the difference in the V-band magnitude be¬ 
tween the two components. The colour and magnitude of 
the secondary component is then simply 


(V - J) b = (V — J) a + A m/a (A3) 

and 


Vb = Va + Am. (A4) 

Hence for TWA 4 we derive, for the two visual compo¬ 
nents A and B, colours and magnitudes of (V — J)a = 
2.399 ± 0.057 mag, Va = 9.442 ± 0.061 mag, (V - J) B = 
2.729 ± 0.057 mag, and Vb = 9.946 ± 0.061 mag. For the un¬ 
certainties on the individual component V — J colours we 
assume the same value as that for the system V — J colour, 
whereas for the individual V-band measurements we add the 
uncertainties on the system V-band and A V-band measure¬ 
ments in quadrature. 


A2 Other cases 

A 2.1 HD 20121 


HD 20121 is a member of Tuc-Hor and consists of a 
triple system (A0V+F7III+F5V according to the Wash¬ 
ington Double Star Catalog; |Mason et al. 2001) in which 
the A+B components are separated by less than 1 arcsec 
and the C component lies a further 3 arcsec away. We 
adopt the combined system V-band magnitude of Vabc = 
5.923 ± 0.005 mag from IMermilliod] ( |2006[ > and a system J- 
band magnitude of Jabc = 5.118 ± 0.030 mag from the 
2MASS PSC, thereby defining a system V — J colour of 
(V — J)abc = 0.805 + 0.030 mag. Tycho-2 provides resolved 
photometry for the combined A+B component and C com¬ 


ponent, from which we derive (using the relation of Ma- 


|majek e~ al. 2006 to transform Vt photometry into John 
son V) a difference between these two visual components 
of AVab.c = 2.642 ± 0.036 mag. Given the early spectral 
types of the components, in conjunction with the age of 
Tuc-Hor, we expect that all three components are on the 


main-sequence. Therefore, to split the photometry we fit an 
empirical (linear) fit to the A/F/G-type main-sequence re¬ 
lation of Pecaut & Mamajek ( 2013[ slope a = 3.311) to de¬ 
rive individual component colours and magnitudes of (V — 
J) AB = 0.758 ± 0.030 mag, (V - J) c = 1.556 ± 0.030 mag, 
Vab = 6.014 ± 0.036 mag, and Vc = 8.656 ± 0.036 mag. 

The reader will note that the spectral types we pro¬ 
vide in Table [l] for the three components of this system dif¬ 
fer from those given in the Washington Double Star Cata¬ 
log. The reason for this is simply that the spectral types 
quoted in the Catalog are inconsistent with the derived 
colours and magnitudes of the individual components. |H 0 g| 
et al. ( 2000[) and Fabricius et al. (2002 I provide resolved 


Tycho-2 (BV )t photometry for the individual components 
of this system. We use the conversion of [Mamajek et al. 


12006) and the trigonometric parallax measurement of 


Leeuwen (2007) to calculate individual absolute My magni¬ 


tudes and Johnson 73 —V colours, and after comparing these 


to the main-sequence relation of Pecaut & Mamajek (20131, 


find that both metrics are consistent with spectral types of 
F4V+F9V+G9V. 


A 2.2 GJ 3322 

GJ 3322 is a member of the BPMG and is described by 


Riedel et al.l (2014) as a triple system (unresolved spec¬ 


tral type of M4IVe) in which the A+C components form 
an unresolved binary and the B component is ~ 1 arcsec 
away. Adopting the combined system V-band magnitude 
of Vabc = 11.504 ± 0.030 mag from UCAC4 and J-band 
magnitude of Jabc = 7.212 ± 0.023 mag from the 2MASS 
PSC, we derive a system V — J colour of (V — J)abc = 


4.292 ± 0.038mag. Riedel et al. (2014) quote a difference in 


magnitude between the combined A+C component and the 
B component of AA'ac.b = 1.03 mag, where the A'-band 
used is that of the CIT system. To estimate what the cor¬ 
responding difference between the two components in the 
V-band is, we first assume that A7\cit — A K s and then 
fit an empirical (linear) fit to the BPMG members in the 
My, V — K s CMD (for which we find a slope of a = 1.665) to 
calculate AVac.b = 2.579 mag (for which we assume an un¬ 
certainty of 0.05 mag). To split the photometry we fit an em¬ 
pirical (linear) fit to the BPMG members (slope a = 1.929) 
to derive individual component colours and magnitudes of 
(V-J) ac = 4.225+0.038 mag, (V-J) B = 5.562+0.038 mag, 
Vac = 11.601 ±0.058 mag, and V B = 14.180 ± 0.058 mag. 


The recent study of Riedel et al. (20141 also includes de- 


blended photometric measurements for the individual com¬ 
ponents, yielding Va = Vc = 12.46 mag, Vb = 13.56 mag, 
(V —J)a = (V —J)c = 4.16 mag, and (V — J)b = 4.91 mag. 

| Riedel et al.| assume that the A and C components are of 
equal-mass, which is clearly not true as |Delfosse et al.| ( |1999| ) 
provides radial velocity amplitudes. On the other hand, our 
method assumes that we are disentangling unresolved bi¬ 
naries, and therefore if one of the components is itself an 
unresolved binary, then our method is insufficiently robust 
to deal with the fact that one of the sources is unnaturally 
brighter because it is a multiple. As we do not have any in¬ 
formation regarding the photometric difference between the 
A and C components, we opt to simply adopt the aforemen¬ 
tioned values of Riedel et ah](|2014j) for the three individual 
components of the system. 
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APPENDIX B: BOOT-STRAPPED COLOURS 
AND QUESTIONABLE MEMBERSHIPS 

B1 Boot-strapped colours 

In Section [272] we discussed cases in which we were forced to 
infer V — J colours for member stars from either the V — K s 
or B — V colour. We used the V — K s colour when only 
the associated 2MASS PSC J-band Qflg was not ‘A’. To 
calculate the uncertainty on the inferred V — J colour, we 
adopted the 2MASS PSC A s -band uncertainty and added 
this in quadrature with the U-band uncertainty. When the 
associated Qflgs of all three 2MASS PSC bandpasses were 
not ‘A’ we used the B — V colour to estimate V — J. In such 
cases, given that the 2MASS PSC photometry is unreliable, 
to assign an uncertainty on the V — J colour, we assumed 
an uncertainty of 0.03 mag combined in quadrature with the 
U-band uncertainty. 

The following stars have V — J colours inferred from 
their V — K B colours: 32 Ori, AK Pic, HD 31647A, 
HIP 88726, HR 136, HR 789, HR 1190, HR 1189, HR 1474, 
HR 1621, HR 2466, HR 6070, HR 7329, HR 7736, HR 8352, 
HR 8911, HR 9016, and HR 9062. 

The following stars have V — J colours inferred from 
their B — V colours: HD 31647B, HR 126, HR 127, HR 674, 
HR 806, HR 838, HR 1249, HR 2020, HR 4023, HR 4534, 
HR 7012, HR 7348, HR 7590, HR 7790, and HR 8425. 


B2 Questionable membership 

In this Section we discuss those stars for which membership 
to a given moving group is either questionable or unlikely. 


available Li diagnostic to infer its youth, we discard GJ 9303 
from our age analysis. 

HR 7329 was classified as a BPMG member by |Zucker^| 
man et al. (20011, however Malo et ah] (20131 demonstrates 
that the radial velocity is discrepant (Av ~ 14kms“ ) with 
regards to that predicted for the BPMG, and that if it is in¬ 
cluded in the prior the probability P v + n becomes zero. |Malo| 


et al. (20131 suggest that either the close-by low-mass com¬ 


panion may affect the systematic velocity or that the radial 
velocity is erroneous due to the fast rotation of the primary 
(usini = 330kms _1 ). With no observations currently avail¬ 
able to categorically rule out its association with the BPMG, 
we retain HR 7329 in the age analysis. 


B2.3 Carina 


AB Pic (HD 44627) has an ambiguous membership having 
previously been assigned to Tuc-Hor bjdZuckerman&Songj 
(20041 only to be revised to Carina b y | Torres et al. (20081. 
The Bayesian analysis of Malo et al.| ( |2013| ) suggests that it 
is more likely a member of the latter (P v + n = 71 per cent 
for Carina compared with only 29 per cent for Columba 
and 0 per cent for Tuc-Hor). The higher probability that it 
belongs to Carina is a result of its Galactic positions YZ 
which are more akin to that of Carina than Tuc-Hor. Fur¬ 
thermore, the reason it has a higher probability of being 
a Columba member (as compared to Tuc-Hor) is because 
there is a difference of only lkms -1 between the predicted 
radial velocities for Columba and Carina (Malo et al.|2013 |. 
Until additional radial velocity measurements are made to 
categorically demonstrate that it is not a member of Carina, 
we retain AB Pic as a member for our age analysis. 


B2.1 AB Dor moving group 

HR 7214 was designated a member of the AB Dor moving 
group by Zuckerman et al. (20111, however Malo et al. (20131 


derives a membership probability (including radial velocity 
and parallax information) of P v +n = 80 per cent for this 
object due to its seemingly anomalous space motion U in 
comparison to the bulk motion of the group as defined by 
bonafide members. We retain HR 7214 for our age analysis 
on the basis that there is no categoric evidence that is should 
be discarded. 


B2.2 BPMG 

HR 789 is a tight spectroscopic binary which also has a 
more distance (25 arcsec) co-moving M dwarf companion 
(2MASS J02394829-4253049). This system was suggested as 
a member of Columba by Zuckerman et al. (2011), however 


Malo et al. (20131 find it is more likely a BPMG member 


(P v + n = 96 per cent). We therefore include HR 789 and 
2MASS J02394829-4253049 as BPMG members for our age 
analysis. 

GJ 9303 (HIP 47133) was proposed as a BPMG member 
by|Schlieder et aT7| ( [20 12a ), however the Bayesian analysis of 
Malo et al. ( 2013| ) suggests that it is more likely associated 
with the field (P„+ 7r = 99 per cent). The main reason for 
this assignment is that its Galactic position in Z is far (~ 
40 pc) from the centre of bonafide BPMG members. With no 


B2.4 Columba 

HD 15115 was previously believed to be a member of the 


BPMG (see e.g. Moor et al. 20061, however Malo et al. (20131 


suggests that it is more likely to be a member of Columba 
(P v +tt = 87 per cent) and therefore we retain this as a mem¬ 
ber of Columba for the age determination. Similarly, whilst 
HD 23524 was originally designated a Tuc-Hor member by 


Zuckerman et al. (2011), Malo et al. (20131 derive a proba¬ 


bility of it being a Columba member of P v + n = 98 per cent, 
and thus we include it in our age determination of Columba. 
In addition, we also include V1358 Ori (BD-03 1386) and 
DK Leo (GJ 2079; originally assigned memberships in Tuc- 
Hor and the BPMG respectively) in our age analysis of 
Columba, as Malo et al. ( |2013 I derive membership prob¬ 
abilities of P v +tt = 85 and 93 per cent respectively. 

The membership status of AS Col (HD 35114) is am¬ 
biguous primarily as a result of a poorly constrained ra¬ 
dial velocity (see e.g. Bobylev, Goncharov, fe Bajkova|2006 


and |Bobylev fe Bajkova 20071 who derive values of 15.2 ± 


1.6kms 1 and 23.9 ± 2.2kms 1 respectively). The analy¬ 


sis of Malo et al. (2013) finds that AS Col belongs to ei¬ 


ther Columba or Tuc-Hor (P v +t = 99 per cent in both 
cases) depending on which radial velocity is adopted. Us¬ 
ing the Bayesian Analysis for Nearby Young AssociatioNs 


(BANYAN; see Malo et al. 2013) for AS Col, but remov¬ 
ing the radial velocity measurement as a prior, we find that 
the star appears to be a bona fide member of Columba. We 
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therefore retain AS Col as a Columba member for our age 
analysis of the association. 

k And (HR 8976) is designated a bona fide member of 
Columba by [Malo et al. (2013) with a probability of P v + n = 
95 per cent. Recent evidence, however, suggests that k And 
is in fact much older (> 100 Myr depending on assumed 
composition; see e.g. Hinklcy ct al. 2013. Brandt & Huang 


20151. Given this evidence we exclude k And from the age 


analysis of Columba. 


B2.5 TWA 


We discard three stars (TWA 14, TWA 18 and TWA 31) 
from the age analysis of the TWA. The mean distance es¬ 
timate of the group is ~ 57 pc with a 1 <t scatter of roughly 


11 — 12pc (Mamajek 20051, however both TWA 14 and 


TWA 18 have distances (trigonometric and kinematic re¬ 
spectively) of > 100 pc and thus if they were members, would 
represent 3<r outliers. In addition, we do not include TWA 31 
in our age analysis due to its designation as a non-member 


by Ducourant et al. (2014). 


B2.6 Tuc-Hor 


GJ 3054 (HIP 3556), HD 12894, HD 200798, and BS Ind 
(HD 202947) are all designated members of Tuc-Hor by 
|Zuckerman fe Song| ( |2004| ), however |Malo et al.| ( |2013| ) 
demonstrated that if their radial velocities are included as 
priors, their respective probabilities of belonging to Tuc-Hor 
significantly decrease due to a difference of ~ 7—8 km s~ x be¬ 
tween the measured and predicted radial velocities. Without 
additional, higher precision, radial velocity measurements 
for these stars, we retain all four in our age determination 
of Tuc-Hor. 

HR 943 was classified as a Tuc-Hor member by |Zuck-| 
|erman et al.| ( |2011[ ) , however |Malo et al.| ( |2013| ) find that 
it is equally likely to be a member of Tuc-Hor or Columba 
(P v + n = 50 and 49 per cent respectively). Until further mea¬ 
surements are able to discriminate between which of the 
moving groups it belongs to, we HR 943 as a member of 
Tuc-Hor for our age analysis. 

HR 6351 and V857 Ara (HD 155915) are given as Tuc- 
Hor members by Zuckerman et al. (2011), and although ex¬ 
hibiting signs of youth (e.g. circumstellar material and Li 


absorption), the analysis of Malo et al. (20131 suggests that 
both are field objects. Without categoric evidence that these 
stars are either field dwarfs or young interlopers belonging 
to another co-moving group, we retain both stars as Tuc-Hor 
members in our age analysis. 


B2.7 Duplicate members 

Given that we have collated memberships from numer¬ 
ous different literature sources, it is possible that stars 
which have been assigned membership to one moving 
group in a particular study, may be classified as a mem¬ 
ber of a different group by another study. We identi¬ 
fied a total of 6 duplicate stars, 5 of which are classi¬ 
fied as members of Tuc-Hor and Columba by |Kraus et al.| 
(|2014|) and |Malo et al.| (|2014a|) respectively, namely CD- 
44 753 (2MASS J02303239-4342232), 2MASS J03050976- 


3725058, 2MASS J04240094-5512223, 2MASS J04515303- 
4647309 and 2MASS J05111098-4903597. In each case the 


star has been assigned membership solely on the basis of 
the measured radial velocity, and for 4 out of the 5 cases 
(excepting CD-44 753) the two independent radial velocity 
measurements of Kraus et al. (20141 and Malo et ah] (2014a I 
agree to within the quoted uncertainties. Given the addi¬ 
tional prior information included in the Bayesian analysis of 
|Malo et al.| we retain these stars as members of Columba, 
but note that the inclusion/exclusion of these stars for either 
group has a negligible effect on the best-fit age. 

TYC 5853-1318-1 (2MASS ,101071194-1935359) was 
also found to be a duplicate, first suggested as a BPMG 
member by Kiss et al. (20111, but more recently advo¬ 
cated as a Tuc-Hor member by Kraus et al. (20141. The 
Bayesian analysis of |Malo et al. ( |2013| ) is unable to cat¬ 
egorically assign membership to only one of the moving 
groups, instead suggesting it could belong to any of the 
BPMG, Tuc-Hor or Columba. Interestingly, [Kraus et al.[ 
(2014) only assigns membership to Tuc-Hor on the basis 
of strong Li absorption i.e. it is young. Its measured ra¬ 
dial velocity (v = 9.3 ± 0.5kms _1 ) is somewhat discrepant 
with respect to that of the bulk of the other Tuc-Hor mem¬ 
bers (Aw = 8.25 kms -1 ), and on the basis of this should 
be considered a non-member (as state d in |Kraus et al.| 
20141. The measured radial velocity of Kiss et al. J201^ 


v = 11.5 ± 1.4 km s ) is consistent with that of 


Kraus 


et al. (20141 to within the uncertainties. Furthermore, both 


of these velocities are consistent with that predicted by|Malo| 


et al. (2013) if it is in fact a member of the BPMG. Based on 


this we prefer to assign membership of TYC 5853-1318-1 to 
the BPMG and include it in our age analysis of the group. 


APPENDIX C: UPDATED r 2 MODEL FOR 
DEALING WITH NON-MEMBER 
CONTAMINATION 

Cl Fitting datasets with non-members 

Cl.l Background 

Non-members whose positions in the CMD lie outside the 
area of the cluster sequence can have an overwhelming effect 
on the fitted parameters, distorting them far away from the 
true values, or make the fit fail entirely. The reason is that 
the total likelihood (the logarithm of which is proportional 
to t 2 ) is the product of the likelihoods that the individual 
data points originate from the cluster sequence. Thus one 
datapoint with a likelihood close to zero for an otherwise 
good fit will drag the model towards it, or even give a prob¬ 
ability of zero for the whole fit. 

Conceptually the most straightforward solution is to 
have a model of the non-nrember contamination in CMD 
space. Perhaps surprisingly it turns out that a very crude 
model is effective; in Bell et al.] |2013| we used a uniform 
distribution over the area delineated by the maximum and 
minimum colours and magnitudes in the dataset (see also 
|van Dyk et al.||2009[ ). In the same paper we showed how 
the uniformly distributed non-member model was formally 
equivalent to a soft-clipping scheme, and in fact actually 
used the latter in the fitting. However, here we will use the 
conceptual framework of a uniform non-member distribu- 
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tion. In part because it is better suited to the problem in 
hand, but also because it allows us to calculate a goodness- 
of-fit parameter, a step which was missing from our earlier 
soft-clipping technique. 

Finally, before embarking on the formalism of this 
method, we should remark that whilst we will talk about 
non-members, in fact we should really refer to stars which 
do not fit our cluster model. For example, equal-mass triples 
will be very slightly above our equal-mass binary sequences, 
and hence in a region of zero probability because our model 
does not include higher order multiples. Therefore, such ob¬ 
jects will be treated by our method as non-members, when 
in fact it is the model which is at fault. In this sense our 
uniform distribution is in part a Jaynes’ fire extinguisher, 
a hypothesis which remains in abeyance unless needed by 


Cl.2 Formalism and implementation 

As in |Naylor fe Jeffrie's] |2006 ) we define the function to be 
minimised as 


data which have a low probability of originating from our 


main hypothesis (see Section 4.4.1 of Jaynes & Bretthorst 


20031. 


= -2 E ln 

i= 1, N 


Ui(c—Ci,m —mi) p(c,m) dcdm (Cl) 


(see also the elegant proof in Walmswell et al. 2013). The 
model of the expected density of stars in the CMD is 
given by p(c, m) and Ui represents the uncertainties for the 
data point i. In terms of Bayes’ theorem, the integral is 
P(M ) P(D\M), and so if we have two competing models this 
should be replaced by P(M n ) P(D\M n ) + P(M C ) P(D\M C ) 
where the subscripts c and n refer to cluster members and 
non-members respectively. If the probability that any given 
star is a member is given by Fi, then P(M n ) = 1 — P; and 
P(M C ) = Fi. Furthermore, P{D\M U ) is the integral of Uip n 
and P{D\M C ) the integral of Uip c , where p c is our usual clus¬ 
ter model and p n is the model of the non-members. If the 
area delineated by the maximum and minimum colours and 
magnitudes in the dataset is A , then where p n is non-zero 
p n = 1/A since it must integrate to one (see |Naylor||2009 ). 
Hence 


t 2 = -2 ^ ln 
i= 1, N 

and thus 


(1 — Fi)p n JJ Uidcdm + Fi JJ Uip c dcdm 

(C2) 


2 


T 


2 ^ In 

i=l, N 


1 -Fi 
A 



Uipcdcdm 


(C3) 


where we have used the facts that where p n is non-zero in 
the CMD it is constant, and that Ui integrates to one. Using 
this formula directly (rather than adding a constant to p c ) 
is very straightforward to implement, since one simply cal¬ 
culates the likelihood in the normal way, adjusts it for the 
probability of membership and adds a constant. This has the 
advantage over our soft-clipping procedure that individual 
stars can be given different membership probabilities. 

A further advantage of this formalism is that we can 


calculate how the position of a star in the CMD modifies 
our estimate of how likely it is to be a member. If we apply 
Bayes’ theorem to the hypothesis M c , that a star is a mem¬ 
ber of the cluster, which we are testing against a dataset D 
then 


P(M C \D) 


P(M C )P{D\M C ) 

P(D) 

P(M C )P(D\M C ) 

P(Mc)P(D\M c ) + P(M n )P(D\M n )' 


(C4) 


Using the same expressions for the probabilities and likeli¬ 
hoods we used for Eqn. |C3] we obtain 


P(M C \D) 


Ft ff Uipcdcdm 

^ + Fi JJ Uidcdm’ 


(C5) 


where in the denominator we have used the same simplifica¬ 
tions as between Eqns. |C2| and |C3| The intuitive interpre¬ 
tation of this equation is that at any point in the CMD the 
ratio of the model densities for members and non-members 
gives the membership probability for a star at that position, 
were its uncertainties in colour and magnitude infinitely 
small. Allowance for uncertainties is made by convolving the 
densities with the uncertainty function. 


C2 Testing the prior membership probabilities 

If the sum of the prior and posterior membership probabil¬ 
ities are very different, this could be an indication that the 
priors were incorrect. A second way of exploring this is de¬ 
scribed in Section 14.21 where we examined the distribution 
of the individual values of r 2 . The cumulative plots in Fig. [ 3 ] 
show how mis-classified members have values of r 2 far ex¬ 
ceeding those predicted by the fitting process. If we lower 
the maximum prior membership probability, then both the 
predicted and measured values of r 2 develop a pedestal at 
the value of r 2 corresponding to the maximum membership 
probability (again shown in Fig. [3|. The simplest way to 
achieve this is to multiply all the priors by a factor, which 
we call P max as it then corresponds to the maximum mem¬ 
bership probability any star can have. We found that ad¬ 
justing P max until either the pedestals in Fig. [3] matched, 
or the sums of the prior and posterior membership proba¬ 
bilities were roughly in agreement gave very similar answers 
for the best value of P max . Given the number of fits we had 
to perform, we took the latter course as straightforward to 
implement in an automated procedure. 

We emphasise that P max should not be adjusted to ob¬ 
tain a reasonable value for the goodness-of-fit [Pr(r 2 )], as it 
is possible have a reasonable value of Pr(r 2 ) but still have 
systematic residuals. This can be seen by examining the cu¬ 
mulative distribution plots for the individual values of r 2 , 
which seem to be a good space in which to examine the qual¬ 
ity of a fit, rather like the residual plots in more conventional 
fitting. Hence there are two metrics for a good fit, the nu¬ 
merical one of the value of Pr(r 2 ), and the qualitative one 
of the match between the predicted and model distributions 
of t 2 . It is possible to adjust P ma x to achieve a good value 
of Pr(r 2 ) , but if this does not achieve a change in shape 
of the r 2 distribution that brings data and prediction into 
agreement, the fit is probably incorrect. 
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C3 Calculation of Pr(r 2 ) 

) showed how to assess whether a 
fit to the data [Pr(r 2 )] by calculat¬ 
ing the probability that a random dataset drawn from the 
model, when fitted to the model would have a value of r 2 
which exceeded that for the observed dataset. The method 
presented for calculating this involved convolving the dis¬ 
tributions of t 2 for each individual datapoint with the dis¬ 
tributions for all other data points. This procedure has the 
advantage of being a numerical equivalent to the way the 
expression for Pr(y 2 ) is derived, but for the data presented 
here there are two significant problems. First the method 
was unacceptably slow for some of the data, this is because 
some of the uncertainties are very large leading to slow con¬ 
volutions. Second the method has to be modified to allow 
for non-members. 


Naylor & Jeffries (2006 


given model was a good 


C3.1 The new method 

We therefore elected to use a more direct method to calcu¬ 
late Pr(r 2 ) where we simply simulated a thousand observed 
datasets using the parameters of the best-fitting model, and 
calculated r 2 for each of them. For the simulation of cluster 
members, the simulated stars should in principle be drawn 
randomly from the stellar mass function multiplied by any 
selection effects, such as magnitude limits. In practice the se¬ 
lection effects are often poorly understood, and so the func¬ 
tion from which the stars are drawn is best defined from the 
magnitude distribution of the real dataset. This problem is 
not unique to this method of calculating Pr(r 2 ), an analo¬ 
gous problem exists for the method described in [Naylor fc] 
|Jeffries| ( |2006| ). 

To calculate Pr(r 2 ) we began by simulating a cluster of 
a million stars with the parameters of the best-fitting model 
(practically one can use the CMD of p created for the fit¬ 
ting process) and grouped the stars into a number of mag¬ 
nitude bins equal to the number of observed data points. 
The boundaries between these bins were set at the mid¬ 
magnitudes between each observed datapoint. We then sim¬ 
ulated a data point from each magnitude bin by first using 
the prior probability of membership to assign it to either the 
cluster or the field. Non-members were simulated by plac¬ 
ing the star randomly within the rectangular area defined 
by the maximum and minimum colours and magnitudes of 
the real dataset (and thus there is a small chance it may be 
placed within the cluster sequence). Cluster members were 
simulated taking their colours and magnitudes from a simu¬ 
lated star randomly chosen from within the magnitude bin. 
Since the bins are more closely spaced where there are many 
real data points, this means the distribution of data points 
in the simulated dataset broadly follows the magnitude dis¬ 
tribution of real data points. Each simulated star was then 
displaced from its model colour and magnitude using the un¬ 
certainties of the corresponding real data point. Following 
this procedure for each magnitude bin resulted in a simu¬ 
lated observation, for which we could calculate t 2 in the 
normal way. The distribution of 1000 values of r 2 calculated 
in this way could then be used to calculate Pr(r 2 ) for the 
real observation. 

Finally, we must allow for the effect of free parameters 
on the expected values of r 2 , since they will be lower for the 



Figure Cl. The probability of a model exceeding a given value 
of r 2 calculated from ten different simulated observations of 30 
stars from a model cluster of a million stars. 


best-fitting model if there are more free parameters than in 
the case where a single model is fitted. For the reasons dis¬ 
cussed in Naylor (2009) we correct for this by scaling the 
values of r" associated with a particular Pr(r 2 ). We first 
subtract the expectation value of r 2 from each of the val¬ 
ues, multiply them by — n), where n here refers to 

the number of free parameters and N to the number of data 
points, and then add back the the expectation value less 
n. Note the scaling factor was incorrectly stated in |Nay-| 
lor] ([2009), and was also incorrect in earlier versions of the 
code, though the corrections are so small that it does not 
significantly affect earlier results. 


C3.2 The accuracy of Pr(r 2 ) 


In principle the cumulative distribution function for the ex¬ 
pected values of Pr(r 2 ) is a function only of the best-fitting 
model and the way in which we select those stars which are 
observed. If we have two different datasets consisting of dif¬ 
ferent stars in the same cluster, provided those stars were 
selected using the same random process, Pr(r 2 ) would be 
the same for both datasets. Unfortunately, the fact we do 
not know the selection criteria for the stars means that, as 


described above (see Appendix C3.11 we have to rely on the 


observations themselves to give us the distribution in magni¬ 
tudes expected for the sample, and so the value we calculate 
for Pr(r 2 ) depends on our sample. This produces an uncer¬ 
tainly in our value of Pr(r 2 ), which it is important to quan¬ 
tify. We performed an example simulation of 30 stars in a 
10 Myr-old clus ter spread over 4 magnitudes in an My, V—J 
CMD. Fig. Cl shows ten of the resulting distributions of r 2 
which show that whilst the shape and the width of the dis¬ 
tribution remain largely constant, the t 2 associated with a 
given probability varies with an RMS (calculated from 100 
models) of ~ 5. For good fits where Pr(r 2 ) ~ 0.5 this is not 
an issue since the 90 per cent width of the distribution is 30 
in t 2 , and so one would almost always conclude the fit is a 
good one. However, in the low-probability tail the test will 
show the fit is probably not a good one, but cannot give a 
useful answer as to how poor it is. 


























